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Abstract 

Context. Accretion-powered X-ray pulsars exhibit significant variability of the Cyclotron Resonance Scattering Feature (CRSF) 
centroid energy on pulse-to-pulse timescales, and also on much longer timescales. Two types of spectral variability are observed. For 
sources in group 1, the CRSF energy is negatively correlated with the variable source luminosity, and for sources in group 2, the 
C^h, opposite behavior is observed. The physical basis for this bimodal behavior is currently not well understood. 

Aims. We explore the hypothesis that the accretion dynamics in the group 1 sources is dominated by radiation pressure near the stellar 
surface, and that Coulomb interactions decelerate the gas to rest in the group 2 sources. 
-4— > ' Methods. We derive a new expression for the critical luminosity, L crit , such that radiation pressure decelerates the matter to rest in 

sources with X-ray luminosity Lx > L cl j t . The formula for L cr i t is based on a simple physical model for the structure of the accretion 
column in luminous X-ray pulsars that takes into account radiative deceleration, the energy dependence of the cyclotron cross section, 
the thermodynamics of the accreting gas, the dipole structure of the pulsar magnetosphere, and the diffusive escape of radiation 
through the column walls. We show that for typical neutron star parameters, L crit = 1.5 x W 31 B 1 ^ 115 erg sec -1 , where B l2 is the surface 
^ . magnetic field strength in units of 10 12 G. 

Results. The formula for the critical luminosity is evaluated for 5 sources, using the maximum value of the CRSF centroid energy 
to estimate the surface magnetic field strength B l2 . The results confirm that the group 1 sources are supercritical (L x > L crit ) and the 
group 2 sources are subcritical (Lx < Lcrit), although the situation is less clear for those highly variable sources that cross over the line 
. L x = L clit . We also explain the variation of the CRSF energy as a consequence of variation of the characteristic emission height with 

luminosity. The sign of this dependence is opposite in the supercritical and subcritical cases, hence creating the observed bimodal 
behavior. 

Conclusions. We have developed a new model for the critical luminosity in accretion-powered X-ray pulsars that explains the bimodal 
dependence of the CRSF centroid energy on the X-ray luminosity Lx. Our model provides a physical basis for the observed variation 
of the CRSF energy as a function of L x for both the group 1 (supercritical) and the group 2 (subcritical) sources as a result of variation 
of the characteristic emission height in the column. 
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1 . Introduction several thousand kilometers out in the accretion disk. The fully- 

ionized accreting plasma is entrained by the magnetic field at the 
X-ray binary pulsars (XRBPs) were first observed by Giacconi et Alfv6n radius> and from there the matter is guided through the 

al. ( 1 97 1 ) and Tananbaum et al. ( 1 972), and now include many mag netosphere, forming accretion columns at one or both of the 
of the brightest sources in the X-ray sky. In XRBPs, the main ma g ne tic poles of the star. As the star spins, the inclination angle 
sequence companion star transfers matter to the neutron star via bet ween the star's magnetic axis and the axis of the accretion 
Roche lobe overflow, or via a strong stellar wind (Frank et al. disk changes, and therefore the Alfven radius varies with the 
2002). The gas forms an accretion disk around the neutron star, S p m p er i od f the star 

and the material spirals inward until the pressure of the star's The observed X -ray emission is powered by the conversion 
dipole magneue field becomes comparable to the ram pressure of gravitationa i pote ntial energy into kinetic energy, which is 

then transferred to the radiation field via electron scattering, and 
ultimately escapes through the walls of the column. The struc- 



of the matter in the disk. This occurs at the Alfven radius, located 
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netic field, with a surface strength B, > 10 12 G, which results in a 
magnetic pressure far exceeding that of either the gas or the radi- 
ation field. The high incident speed of the freely-falling plasma, 
~ 0.6 c, creates very high temperatures, T ~ 10 8 K. However, the 
observed X-ray pulsar spectra are highly nonthermal, indicating 
that the accreting gas is unable to equilibrate during the accretion 
timescale. In this situation, bulk and thermal Comptonization 
play key roles in establishing the shape of the observed spectra 
(Becker & Wolff 2007). 

The X-ray spectra of many XRBPs contain cyclotron reso- 
nant scattering features (CRSFs) appearing as absorption lines. 
The features are caused by resonant scattering of photons off 
plasma electrons whose energy is quantized according to their 
Landau level (see e.g. Triimper et al. 1978; Isenberg et al. 1998; 
Araya-Gochez & Harding 2000). The CRSFs, when detected, 
provide a direct measurement of the magnetic field strength at 
the characteristic altitude of the X-ray emission. The energy of 
the fundamental line and the spacing between the harmonics are 
approximately proportional to the B-field strength. 

Many XRBPs display X-ray spectra that vary significantly 
with luminosity on timescales much longer than the pulsation 
period. In particular, variations in the energy of CRSFs as a func- 
tion of luminosity on timescales of days to months have been 
detected in V 0332+53 (Mowlavi et al. 2006; Tsygankov et al. 
2010), 4U 01 15+63 (Mihara et al. 2004; Tsygankov et al. 2007), 
and HerX-1 (Staubert et al. 2007; Vasco et al. 201 1). In addition 
to the longer-term variability, there is also mounting evidence 
for pulse-to-pulse variability, in which the spectral hardness, the 
centroid energy of the CRSF, and the luminosity vary in a corre- 
lated way (Klochkov et al. 2011). This short-timescale variabil- 
ity is likely related to the non-stationarity of the accretion flow, 
perhaps indicating that the entrainment of matter from the disk 
onto the magnetic field lines results in filaments and blobs of 
accreting gas which are then channeled onto the star in a non- 
uniform way. 

The data from both long-term and short-term (pulse-to- 
pulse) observations point to the existence of two types/modes 
of spectral variability (see discussion by Klochkov et al. 2011). 
For sources in group 1 (e.g., V 0332+53), the centroid energy of 
the CRSF is negatively correlated with luminosity. For sources 
in group 2 (e.g., HerX-1), the opposite behavioris observed. The 
type of spectral variability is likely driven by the mode of accre- 
tion, which in turn is determined by the luminosity (see discus- 
sion in Staubert et al. 2007). Staubert et al. (2007) and Klochkov 
et al. (201 1) have proposed that for sources in group 1, the decel- 
eration of the flow to rest at the stellar surface is accomplished by 
the pressure of the radiation field, and in the group 2 sources the 
deceleration occurs via Coulomb interactions. In this interpre- 
tation, a given source falls in one group or the other depending 
on the value of its X-ray luminosity, Lx, relative to the critical 
luminosity, L cr j t . The hypothesis is that the group 1 sources are 
supercritical, with X-ray luminosity Lx > L cr i t , and the group 2 
sources are subcritical (Lx < Lent)- 

The theory predicts that sources in their supercritical state 
should display a negative correlation between the luminosity and 
the cyclotron energy, while sources in the subcritical state should 
display the reverse behavior. Geometrically, the variation of the 
CRSF energy with luminosity is connected with variation of the 
characteristic emission height, which is the altitude in the accre- 
tion column where the cyclotron absorption feature is imprinted 
on the observed spectrum. The variation of the emission height 
as a function of luminosity in the subcritical and supercritical 
cases is indicated schematically in Fig. 1. 



The general picture described above provides a qualitative 
basis for the interpretation of the observed correlated variation 
of the CRSF centroid energy with X-ray luminosity in some 
XRBPs. However, in order to obtain a quantitative understanding 
of these observations, one must develop a more detailed physical 
model for the critical luminosity, and for the dependence of the 
CRSF energy on the luminosity in the subcritical and supercriti- 
cal regimes. The first goal of this paper is to derive a new expres- 
sion for the critical luminosity, taking into account the magneto- 
spheric connection between the radius of the accretion column 
and the Alfven radius in the disk, and the energy and angle de- 
pendence of the cyclotron scattering cross section. The second 
goal is to examine the dependence of the CRSF centroid energy 
on the luminosity in the subcritical and supercritical sources. 

The remainder of the paper is organized as follows. In 
Sect. 2, we obtain a fundamental expression for the critical lumi- 
nosity that depends on the stellar mass, radius, and surface mag- 
netic field strength. In Sect. 3 we develop simple physical models 
for the variation of the characteristic emission height as a func- 
tion of the luminosity for subcritical and supercritical sources. 
In Sect. 4, we evaluate the critical luminosity based on measure- 
ments of the CRSF centroid energy for several XRBPs. We use 
our subcritical and supercritical models for the variation of the 
emission height to predict the variation of the CRSF energy as a 
function of luminosity. The predicted spectral variability is com- 
pared with the observational data for each source. We discuss our 
results and draw conclusions in Sect. 5. 

2. Critical luminosity 

2.1. Eddington luminosity and radiative deceleration 

First we recall the definition of the standard Eddington luminos- 
ity, Lead, for spherically symmetric accretion onto a central mass 
M*. If the accreting gas is fully-ionized hydrogen, we obtain 

AnGM„m v c 
L Ed d = , (1) 

CTj 

where cr T is the Thomson cross section, m p is the proton mass, 
c is the speed of light, and G is the gravitational constant. When 
the X-ray luminosity Lx = Lem, the rate at which momentum is 
transferred to the gas via Compton scattering balances the gravi- 
tational force. Hence if Lx > Le<j<i> the net force is in the outward 
direction and the gas decelerates as it falls toward the central 
mass. 

We need to make two adjustments to Equation (Q]) in order 
to compute the effective Eddington limit, L* dd , appropriate for 
treating X-ray pulsar accretion columns. The first adjustment is 
to replace the Thomson cross section cr x with cry, which repre- 
sents the mean scattering cross section for photons propagating 
parallel to the magnetic field axis. The second adjustment is to 
take the geometry of the accretion flow into account by reducing 
the luminosity by the ratio of the column cross-sectional area di- 
vided by the surface area of the star. Employing these corrections 
yields for the effective Eddington limit 

0- T 7ir^ GM,m v c nr^ 

^dd = L Edd — — ~2 = — zj- > ( 2 ) 

cry 4nRi CT|| Rj 

where R t is the stellar radius and ro denotes the radius of the 
accretion column, which we assume to have an approximately 
cylindrical geometry. 

Based on Equation (0, Basko & Sunyaev (1976) concluded 
that for X-ray luminosities Lx i 10 36 erg s _1 , the incident, freely- 
falling gas is decelerated by a vertical flux of radiation that is 



2 



P. A. Becker et al.: Bimodal spectral variation in X-ray pulsars 



locally super-Eddington. The scattering of the incident radiation 
removes kinetic energy from the electrons (and from the protons 
via Coulomb coupling), thereby decelerating the gas. Although 
the radiation flux seen by the gas inside the column is propa- 
gating upward, the X-rays that ultimately carry away the kinetic 
energy actually escape through the walls of the column, rather 
than the top, unless the luminosity Lx 5 10 36 erg s _1 . 

The deceleration of the gas begins when the freely-falling 
material encounters a radiation-dominated shock whose height 
above the star increases with increasing luminosity, reaching an 
altitude of several kilometers for Lx ~ 10 37 ~ 38 erg s~'. As the gas 
passes through the shock, the accretion velocity is reduced by a 
factor of ~ 7. Even though the radiation flux inside the column 
is super-Eddington, the height of the radiation-dominated shock 
is stable. This reflects the fact that the shock is a wave structure, 
and is not composed of a fixed population of particles. Matter 
moves through the shock and decelerates, but the shock height 
remains fixed, unless the luminosity changes. Unlike a classical 
gas-mediated shock, the radiation-dominated shock is not dis- 
continuous, and instead has a thickness that is a few times larger 
than the mean-free path for radiation scattering. 

Below the radiation-dominated shock, the matter is further 
decelerated in the hydrostatic "sinking regime," in which the 
remaining momentum is transferred to the radiation field and 
radiated away through the column walls (Basko & Sunyaev 
1976). The specific mechanism accomplishing the final decel- 
eration to rest at the stellar surface in the sinking region de- 
pends on the luminosity of the accretion flow (see Fig. 1). At the 
highest luminosities, Lx ~ 10 37 ~ 38 erg sec -1 , the radiation field 
accomplishes the deceleration all the way down to the stellar 
surface (Basko & Sunyaev 1976). At intermediate luminosities 
Lx ~ 10 36 ~ 37 erg sec -1 , the final phase of deceleration may occur 
via Coulomb breaking in a plasma cloud just above the stellar 
surface (Nelson et al. 1993). It is expected that at very low lumi- 
nosities, Lx i 10 34 ~ 35 erg sec -1 , there is no radiation-dominated 
shock at all, and the material passes through a conventional gas- 
mediated shock at altitude z = h g before striking the stellar sur- 
face (Langer & Rappaport 1982). 

The angular pattern of the emitted radiation also depends 
on the luminosity (see Fig. 1). In high-luminosity sources 
(Lx ~ 10 37 ~ 38 erg s~'), the emitted radiation primarily es- 
capes through the column walls in the sinking region, form- 
ing a "fan beam" (Davidson 1973). For low-luminosity sources 
CLx f 10 35 erg s _1 ), the emission escapes from the top of the col- 
umn, forming a "pencil beam" (Burnard et al. 1 99 1 ; Nelson et al. 
1993). In the intermediate range, L x S 10 35 ~ 37 erg s _1 , the emis- 
sion pattern may be a hybrid combination of these two types 
(Blum & Kraus 2000). 

Focusing on the high-luminosity case for now, we can esti- 
mate the luminosity required to decelerate the gas to rest at the 
stellar surface by considering the physical processes occurring 
in the sinking region below the radiation-dominated shock. The 
accreting matter approaches the top of the shock with the inci- 
dent free-fall velocity, which we approximate using the value at 
the stellar surface, 



v ff = 



(2GM t 



1/2 



(3) 



Advection is dominant over diffusion in the shock, and therefore 
very little radiation energy escapes through the walls of the ac- 
cretion column in the vicinity of the shock (Bumard et al. 1991). 
Hence the shock jump conditions are well approximated by the 
standard Rankine-Hugoniot relations for a gas with adiabatic in- 
dex y - 4/3 (Basko & Sunyaev 1976). In this case, the matter 



leaves the shock with the post-shock velocity 



1 



1 (2GM, 



Vps = - Vff = - 



7 



7\ R, 



1/2 



(4) 



where we have assumed that the shock is strong, which is rea- 
sonable in the luminous sources (Becker 1998). If the altitude 
of the radiation-dominated shock above the stellar surface is H, 
and the gas decelerates at a constant rate a from the post-shock 
velocity Vff/7 to rest at the stellar surface in the dynamical time 
fdyn, then we can write the simple kinematical relations 



H = \ at lyn 



a ?d 



yn 



(5) 



Upon elimination of f^yn, we obtain for the required upward ac- 
celeration 



'ps 



GM, 



2H 49RM 



(6) 



Since the effective gravity is reduced by the pressure of the 
radiation field, the net acceleration can also be related to the lu- 
minosity Lx via 



Lx 
L 



Edd 



GM, 



Setting Eqs. (O and © equal and solving for Lx yields 



Lent — L c 



Edd \49H 



+ 1 



(7) 



(8) 



Substituting for L* dd using Eq. ©, we obtain for the critical lu- 
minosity 



Lcrit — 



GM,m v c nrf. i R t 



Rl \ 49 H I 



(9) 



Our goal is to express the parameters ro, a\, and H appearing on 
the right-hand side of Eq. (01 in terms of observable quantities. 

2.2. Radiation-dominated shock height 

The altitude, H, of the radiation-dominated shock can be esti- 
mated by considering the relationship between the dynamical 
timescale for deceleration, fd yn , and the photon escape timescale, 
fdiff, which is the mean time it takes the photons to diffuse 
through the walls of the accretion column. In order for the gas 
to come to rest at the stellar surface, these two timescales must 
be comparable in the sinking region below the shock, which is 
a general property of accretion flows onto white dwarf stars and 
neutron stars (e.g., Imamura et al. 1987). Combining Eqs. © 
and (0, we obtain for the dynamical time 



R 



2H 

fdyn = = 14//L 

v ps \ 2GM t 



1/2 



(10) 



The escape timescale for the photons to diffuse through the col- 
umn walls is estimated by writing 



..diff ' 



vf = 



— , t ± = r n e cr T 



(ID 



where n e is the electron number density, is the photon dif- 
fusion velocity perpendicular to the column axis, and t ± is the 
perpendicular optical thickness. The Thomson cross section cr T 
is appropriate for photons propagating perpendicular to the col- 
umn axis (Wang & Frank 1981). The electron number density n e 
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Figure 1. Schematic illustration of the geometry of the accretion column and the variation of the characteristic emission height and 
emission beam pattern with increasing luminosity Lx- (a) subcritical, Lx < L cou \ < L a i t , pencil beam; (b) subcritical, Lx 5 L cou i < 
Lent, pencil beam; (c) subcritical, L cou i < Lx ~ L cr \ t , intermediate beam pattern; (d) supercritical, Lx J L c rit, fan beam. 



appearing in Eq. dTTT > can be eliminated using the mass conser- 
vation relation, 

(12) 



M — 7rr Q n e m p v . 



where M is the accretion rate and v is the inflow velocity, defined 
to be positive. Combining relations, we can express the escape 
time through the walls as 



M(Tt 



(13) 



R*/(49H) <sc 1 for sources close to or above the critical lumi- 
nosity, and therefore Eq. ((9) reduces to 



Lent — 



GM*m p c 7ir ( , 



Ri 



(17) 



in agreement with Burnard et al. (1991). Note that in this limit, 
the critical luminosity simply reduces to the effective Eddington 
value given by Eq. ©. 



The deceleration in the sinking region begins on the down- 
stream side of the shock, and therefore we set v = v ps in Eq. ( TOT l 
and equate fd yn and f esc to obtain 



H ■■ 



Mcr T 
2nm v c 



(14) 



which is essentially the same result obtained by Burnard et al. 
(1991). Expressing the accretion rate in terms of the luminosity 
using the relation 

U - ^ (15) 



R 



yields the equivalent expression 

-l 



H = 1.14x10 s cm 



M„ 



I AM, 



(lOkrjflO^ergsec- 1 ) ' (16) 



This confirms that the shock is located a few kilometers 
above the stellar surface in the luminous sources with Lx ~ 
10 37-38 erg s -i ( Basko & Sunyaev 1976). It follows that 



2.3. Connection between column radius and Alfven radius 

In this section, we wish to relate the critical luminosity L cr ; t in 
Eq. (TTTb to the magnetic field strength at the stellar surface, fi*, 
by utilizing the connection between the radius of the accretion 
column, ro, and the Alfven radius in the disk, R&- The inclina- 
tion angle between the axis of the accretion disk and the star's 
magnetic axis varies with a period equal to the pulsar's spin pe- 
riod, and this causes an associated variation of the Alfven radius. 
However, for our purposes here, an adequate approximation is 
obtained by using Eq. (13) from Lamb etal. (1973), which yields 



R A = 2.73 x 10 7 cm — 



A 



O.1/\1.4M 

B* ^ 



UOkrrJ 



x 



\10 12 G/ 



, -2/7 



10 37 erg sec 



(18) 

where the constant A = 1 for spherical accretion and A < 1 
for disk accretion. A variety of uncertainties are folded into A, 
such as the spin-averaging of Rp, and the possible role of plasma 
shielding and other magnetospheric effects. Based on Eq. (2) 
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from Harding et al. (1984), A can be approximated in the disk 
application using 

A as 0.22 q- i8/69 , (19) 

where a < 1 denotes the Shakura-Sunyaev parameter (Shakura 
& Sunyaev 1973). Although it is difficult to estimate a with any 
certainty, we generally expect to find a ~ 0.01 - 0.1. We there- 
fore set A = 0.1 in our numerical applications. 

The Alfven radius in the disk is connected with the outer 
surface of the accretion column through the dipole shape of the 
pulsar magnetosphere. The equation for the shape of the criti- 
cal field line as a function of the polar angle 6 is given by the 
standard dipole formula 



R = R A sin 2 8 . 



(20) 



Setting the radius R equal to the stellar radius R* yields for the 
critical angle at the outer edge of the accretion column 



A' 
R~a 



(21) 



Using the small-angle relation 6 C » sin# c , we obtain for the col- 
umn radius 

RA (22) 



By substituting for the Alfven radius in Eq. d22l using 
Eq. (fT8l . we find that the expression for the column radius ro 
can be rewritten in cgs units as 



r = 1.93 x 10 5 cm 



ay 112 ! m* r 1/I4 / R. \ 11/14 



0.1/ \1.4M 

B. \ 2n 



UOkm/ 



(io 12 g) 



10 37 erg sec 



1/7 



(23) 



Using Eq. (1231 to substitute for ro in Eq. ( TTTl i. and setting Lx 
L CI it, we obtain for the critical luminosity the new expression 



t \-7/5 / \-7/S 



L crit = 7.79 x 10 35 erg sec" 1 — — 



0.1/ \cr T 



lAMo 



( — — — ) 3 5 ( — — — ) 

1 10km/ \10 12 G/ 



-4/5 



(24) 

where we have also introduced the Thomson cross section cry 
as a convenient scaling for the parallel scattering cross section 
cr\\. The next step is to evaluate the cross section ratio cr||/<x T in 
terms of observable source parameters. 

2.4. Electron scattering cross section for parallel propagation 

In typical X-ray pulsars, most of the observed radiation is emit- 
ted at energies below the cyclotron energy, E cyc . Hence the cross- 
section ratio ct||/ctt can be roughly approximated using (e.g., 
Arons et al. 1987) 

C"T V^cyc/ 

where £ is a measure of the mean energy of the photons prop- 
agating parallel to the magnetic field, and the cyclotron energy 
E cyc is given by 



/, u -U.58keV( |( 



\10 12 G/ 



(26) 



The mean photon energy E in Eq. ( 1251 can be estimated 
observationally by integrating the spectrum for a given source. 
However, in luminous X-ray pulsars, most of the observed radi- 
ation escapes through the walls of the accretion column, perpen- 
dicular to the magnetic field, and therefore the observed spec- 
trum may not be representative of the distribution of photons 
propagating along the column axis. As an alternative, we can 
estimate E based on the thermal structure of the accreting gas. 
Specifically, we assume that 



E = wkT, 



eft' ■ 



(27) 



where T e s is the effective temperature of the radiation in the 
post-shock region, k is Boltzmann's constant, and the constant 
w depends on the shape of the spectrum inside the column. 
We expect that w ~ 1-3, with the lower value correspond- 
ing to bremsstrahlung and the upper value to a Planck spectrum. 
Detailed models suggest that the spectrum inside the column is 
dominated by bremsstrahlung emission (Becker & Wolff 2007), 
and therefore we will set w = 1 in the numerical results pre- 
sented later. 

The effective temperature is related to the post-shock radia- 
tion pressure, P r , via 



oT* = 3P r . 



(28) 



The value of P r can be estimated using the momentum balance 
relation 

, Mvff 

P r = Pff v 2 =^, (29) 
nr 

where ptf and vjj (Eq. (0) denote the upstream mass density and 
velocity, respectively, just above the shock. Eliminating P T be- 
tween Eqs. ( l28l i and (1291 . and substituting for M and ro using 
Eqs. (Q~5]! and (J23j, we find that 



r eff = 4.35 x 10' K 



A \ 1/4 / M. 



0.1/ \1.4M 



-5/56 



UOkm/ 



-15/56 



(io i2 g) 



\5/28 



10 37 erg sec 



(30) 

Combining Eqs. ( l25l ). ( 1261 ). d27l >. and (1301 ). we obtain for the 
required cross section ratio the result 



£7-11 / A 

-- =0.106 — 

O-T I 0.1 



1/2 



M,. 



I AM, 



5/28 / r. y 

llOkm/ 



-15/28 



(io 12 g) 



■12/7 



,5/14 



10 37 erg ; 



(31) 



Using this result to substitute for cr\\l<jj in Eq. (l24l) yields the 
final expression for the critical luminosity as a function of the 
surface magnetic field strength, 



Lent = 1.49 X 10 37 erg sec" 1 (A 



-7/5 



-28/15 



M, 



1.4M 



29/30 



/ R t \ 1,w ( B„ \ 
1 10 km/ \10 12 G/ 



16/15 



(32) 



For typical neutron star parameters, with M, = 1 .4 M , R* = 
10km, A = 0.1, and w = 1, we obtain L crit = 1.49 X 
10 37 erg sec^'fij^ 15 ' where Bn is the surface magnetic field 
strength in units of 10 12 G. In Sect. 4 we plot the critical luminos- 
ity and compare it with the variability data for several XRBPs. 
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3. Variation of emission height 

The new expression for the critical luminosity given by Eq. (1321 
allows us to separate accretion-powered X-ray pulsars into sub- 
critical and supercritical categories. Our hypothesis is that in the 
subcritical sources with variable luminosity Lx, the cyclotron 
energy E cyc will exhibit a positive correlation with Lx, and in 
the supercritical sources the reverse behavior will be observed. 
The observed CRSF is imprinted on the spectrum at the altitude 
where most of the emitted radiation escapes from the accretion 
column. In order to quantify the expected behaviors in the sub- 
critical and supercritical regimes, we must therefore explore the 
expected variation of the emission height as a function of lumi- 
nosity for both types of sources. The geometry of the super- and 
subcritical sources is depicted schematically in Fig. 1. 

3. 1 . Supercritical sources 

In the supercritical sources (luminosity Lx J L C ntX radiation 
pressure dominates the flow dynamics all the way to the stel- 
lar surface. Inside the radiation-dominated shock, the infalling 
matter begins to decelerate as it first encounters the "cushion" 
of radiation hovering at the shock altitude (Davidson 1973). At 
this altitude, there is a local balance between downward advec- 
tion and upward diffusion of radiation, and therefore the photon 
distribution is roughly static. Most of the kinetic energy of the 
accretion flow is carried away by the scattered radiation, which 
is likely to be beamed in the downward direction due to special- 
relativistic aberration (e.g., Ferrigno et al. 2009). Below the 
shock altitude, the photons are trapped by advection, although 
they eventually manage to escape by diffusing through the walls 
of the column. Hence the observed radiation does not escape 
from the shock altitude H, but rather from a lower altitude. 

Our goal here is to estimate the typical altitude, denoted by 
h s , at which the photons diffuse through the column walls to 
form the observed X-ray spectrum in the supercritical case. We 
assume that the observed CRSF is imprinted at this altitude, be- 
cause at higher altitudes the photons have not had enough time 
to diffuse through the column and escape through the walls. 
Conversely, at lower altitudes, the increasing density of the gas 
in the column inhibits the escape of radiation. We therefore ex- 
pect that the CRSF energy will reflect the cyclotron energy at the 
altitude z — h & where the final deceleration phase begins. 

We can estimate the emission height /i s in the supercritical 
sources by ensuring that all of the kinetic energy is radiated 
through the walls by the escaping photons in the altitude range 
< z < h s (Basko & Sunyaev 1976). Working in the frame co- 
moving with the plasma in the accretion column, we note that 
the fraction of the radiation escaping through the walls in the 
comoving time interval dt' is equal to dt' /f esc , where f e sc is the 
escape time given by Eq. ( fT3b . The requirement that all of the ra- 
diation escapes by the time the matter reaches the stellar surface 
is therefore expressed by the integral condition 



Jo I dz' 



dz' 

fesc(z') 



= 1 . 



)o I dz' 

Using Eq. dT3l > to substitute for f esc yields 
' hs nm p vc dz' 



MlTj V e ff 



1 • 



(33) 



(34) 



where the effective velocity for the photon transport is defined 
by 

dz' 



Veff = 



dt' 



The flow is expected to be almost perfectly "trapped" in 
the region below z' = h s , meaning that advection and diffu- 
sion are nearly balanced, leaving very little net flux of radia- 
tion (Becker 1998). This implies that the effective velocity v e ff is 
much smaller than the flow velocity v. We define the parameter 
£ as the ratio of these two velocities, 



_ Veff 
V 



(36) 



We demonstrate in Appendix A that the value of £ can be esti- 
mated using 

^ = 77^' (37) 

where Aloo denotes the incident (upstream) Mach number of the 
flow with respect to the radiation sound speed. In the sinking 
region below the shock, the effective velocity approaches zero as 
the gas settles onto the stellar surface. The relationship between 
the upstream Mach number A1 M and the X-ray luminosity Lx is 
plotted in Fig. 12 from Becker (1998). For the parameter range 
of interest here, it is sufficient to adopt a constant value for £ 
in the range £ ~ 10 -2 - 10 -3 . The low value for the effective 
velocity tends to make the emission region more compact in the 
supercritical sources. 

Combining relations, we find that 



f 

Jo 



■dz' = 1, 



lo Mo- T % 

and therefore the altitude of the emission region is given by 



In 



mn p c mn p cGM t 



(38) 



(39) 



where the final result follows from Eq. (TT5t . We can also express 
h s in cgs units using 



h s = 2.28 x 10 3 cm (—) I M * ) 
\0.017 \1AM B ) 

X (l0km)( 10 37 erg sec 1 / 



(40) 



Note that the emission height in the supercritical sources varies 
in proportion to the luminosity Lx (see Fig. 1). 
Based on Eqs. dT~6b and d40b . we conclude that 



^=2.0£«1 



(41) 



(35) 



and therefore the characteristic height of emission in the super- 
critical sources is located far below the altitude of the radiation- 
dominated shock. Eq. ( l40b indicates that the height of the emis- 
sion region h s scales in proportion to the luminosity Lx in the 
supercritical sources, which is consistent with the observed be- 
havior in the group 1 sources (Klochkov et al. 201 1). 

3.2. Subcritical sources 

In the subcritical sources (luminosity Lx S L c ritX the matter 
still passes through a radiation-dominated shock, which accom- 
plishes the initial deceleration, but the pressure of the radiation 
is insufficient to bring the matter to rest at the stellar surface 
(Basko & Sunyaev 1976). In this case, the final stopping occurs 
via direct Coulomb interactions close to the base of the accre- 
tion column (Burnard et al. 1991). Our goal in this section is 
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to estimate the characteristic emission height for the subcritical 
sources, denoted by h c , which is the altitude at which Coulomb 
interactions begin to decelerate the plasma to rest. The emis- 
sion is expected to be concentrated in this region because this 
is essentially the first opportunity that the radiation inside the 
column has to diffuse through the walls. At lower altitudes, the 
radiation is trapped in the column due to the increasing density, 
and therefore the emission through the walls tapers off as z — » 0. 
Hence we expect that in the subcritical case, the CRSF energy 
is imprinted at the altitude z = h c where the strong Coulomb 
deceleration begins. 

The Thomson optical depth, t», required to stop the flow via 
Coulomb interactions can be estimated in the typical pulsar mag- 
netic field regime using Eq. (3.34) from Nelson et al. (1993) to 
write 

M, \ 2 /^)- 2 ^_ ; 
1 10 km/ ' " 



51.4 



I AM, 



ln(2n max ) 



where the maximum excited Landau level, n max , is given by 



2E, 



(43) 



eye 



A summary of the derivation leading up to Eq. (l42l is provided 
in Appendix B. Adopting typical X-ray pulsar parameters, we 
find that r* ~ 20, which is the value utilized in our numerical 
examples. 

We can use Eq. ( l42l to estimate the emission height in the 
subcritical sources, /i c , as follows. The Thomson depth r as a 
function of the altitude z measured from the stellar surface is 
computed using 



T(Z) 



F p(z')o-t 
Jo m v 



dz' 



(44) 



where p = n e m p is the mass density, given by (see Eq. (fT2l ) 



P(z) = 



M 



nr^viz) 



(45) 



Assuming that the gas decelerates uniformly in the Coulomb 
stopping region (starting at altitude h c ) from the post-shock ve- 
locity v ps , we find that the required deceleration is given by 
a - vj; s /(2/i c ) (cf. Eq. (0). The velocity profile v(z) associated 
with the constant deceleration a is computed using the standard 
kinematical relations 



dv dv 1 dv 
dt dz 2 dz 



(46) 



where the negative sign appears because we have defined v and a 
to be positive. Setting v ps = Vff/7 (see Eq. (0]l), we obtain for the 
deceleration a = v^/(98/i c ). Upon integration of Eq. d46i l. we 
therefore find that the velocity profile in the Coulomb stopping 
region is given by 



v(z) = — J — , 
7 V h c . 



(47) 



where vs is evaluated using Eq. (01. 

Using Eq. ( |47| ) to substitute for v(z) in Eq. d45l > and carrying 
out the integration in Eq. (l44l . we obtain for the optical depth 
profile 

UMcr T (2GM„Y 1/2 , — 
t{z) = — 5— !- — — yfhTz. (48) 



Finally, setting z — h c and t = t», we find that the Thomson 
optical depth required for Coulomb stopping is given by 



T* = 



14Mcr T /i c / 2GM, 
ntinu \ R* 



-1/2 



(49) 



which can be rearranged to obtain for the subcritical emission 
height 

UMo-t \ R* ) 

Substituting for M and ro using Eqs. < fT3T > and d23l . respectively, 
yields the equivalent cgs expression 

19/14 



\0.lJ \2O/\1.4M / VIOkm/ 

x (io^g) 



10 37 erg sec 

(51) 

This result indicates that the emission height in the subcritical 
sources decreases with increasing luminosity, which is consis- 
tent with the behavior observed in the group 2 sources (Staubert 
et al. 2007). As indicated in Fig. 1, in the subcritical case, a de- 
crease in the luminosity causes the beam pattern to transition 
from a pure fan configuration to a hybrid pattern that includes 
a pencil component. At very low luminosities, the shock essen- 
tially sits on the stellar surface, and the emission occurs via the 
pencil component only (Burnard et al. 1991; Nelson et al. 1993). 

3.3. Transition from Coulomb stopping to gas shock 

Eq. dBll l gives the height of the emission zone in the subcriti- 
cal case, based on the assumption that the final deceleration to 
rest at the stellar surface occurs via Coulomb interactions. We 
can estimate the minimum luminosity (or accretion rate) such 
that Coulomb interactions are capable of stopping the flow by 
requiring that h c < H, where H is the altitude of the radiation- 
dominated shock. If this condition is violated, then the incident 
velocity of the gas entering the Coulomb deceleration region be- 
comes essentially the full free-fall velocity, Vff, rather than the 
reduced post-shock velocity, v ps = vs/1 . This in turn causes the 
Thomson depth r to fall below the value t» ~ 20 required for the 
gas to be effectively stopped via Coulomb interactions. It is not 
completely clear what happens in this case, but we expect that 
the final phase of deceleration will occur via passage through a 
gas-mediated shock near the stellar surface (Langer & Rappaport 
1982). 

By combining Eqs. (fT~4T >. (TT~5b and 1501 . we can show that the 
condition h c < H implies that 

Substituting for ro using Eq. (1231 and rearranging, we obtain 



Lx > L cou \ 



(53) 



where 



L OOU1 = 1.17X10- erg sec- (A) 7/12 (lip (_^_) 

/ R, Y lV2A l B t \-'/ 3 
X 1 10km/ U0 12 G/ 



11/8 



(54) 
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For luminosities Lx ~ Lcouh we expect that the characteristic 
emission height settles down onto the stellar surface. For very 
low luminosities, Lx ~ 10 34 ~ 35 erg sec -1 , the radiation-dominated 
shock and the Coulomb atmosphere both dissipate, and the 
matter strikes the stellar surface after passing through a gas- 
mediated shock (Langer & Rappaport 1982), as indicated in 
Fig. 1. 



4. Applications 

Our final result for the critical luminosity as a function of the 
surface magnetic field strength B* is given by Eq. (l32l . which 
can be rewritten as 



Led, = 1.28 x 10 37 erg sec 



A 



-7/5 



W 



-28/15 



M. \ 29/30 / R* \ 1/w / E. v 16/15 



I AM, 



1 10 km/ UOkeV/ 



One can see that the sign of the correlation between the cy- 
clotron line energy and the luminosity, when detected, is con- 
sistent between the long-term and pulse-to-pulse studies. This 
suggests that both are reflecting the same underlying physics. It 
is clear that the variation of L cyc with Lx for sources on each 
side of the critical line Lx = L cr j t is qualitatively consistent 
with the theoretical predictions, i.e., the correlation between L cyc 
and Lx is negative for the supercritical sources, and positive for 
the subcritical ones, reflecting the positive variation of the emis- 
sion height h = h s (Eq. ( 1401 ) with increasing Lx for the former 
sources and the negative variation of the emission height h = h c 
(Eq. ( BTT i) for the latter sources. In the next section, we explore 
the variation of L cyc as a function of L x in more detail. 

4. 1 . Variation of L cyc with L x 

(55) The observed value of £ cyc is connected with the local field 
strength B at the emission altitude h via Eq. ( f26l >. where B has 
the dipole dependence 



where 



11.58keV 



(io 12 g) 



(56) 



denotes the surface value of the CRSF energy. This relation is 
indicated by the dashed blue line in Fig. 2. Sources to the right 
of this line are radiating supercritically, and consequently radia- 
tion pressure accomplishes the deceleration all the way down to 
the stellar surface. For sources to the left of this line, the final 
deceleration occurs via Coulomb interactions. 

It is now interesting to compute the critical luminosity for 
a number of sources and to compare it with the observed vari- 
ability of £ cy c as a function of Lx for subcritical and supercriti- 
cal cases. We have selected for this study the sources for which 
the behavior of the cyclotron line energy with luminosity has 
been studied sufficiently well, including both the variations on 
long timescales (days to months) and short timescales ("pulse- 
to-pulse"). Our source sample is listed in Table 1, where we 
also summarize the reported L cyc -Lx behavior (positive/negative 
correlation or no dependence) for each source. The two panels 
in Fig. 2 depict the variability of E cyc with luminosity both on 
"pulse-to-pulse" (left) and longer (right) timescales for all of the 
sources considered here. The corresponding references for the 
data are given in Table 1 . The data for the longer-term variability 
are based on studies of the phase-averaged emission conducted 
either during outbursts of the transient sources (Tsygankov et al. 
2007, 2010), or following the long-term variation of the emis- 
sion from the persistent sources (Staubert et al. 2007). For the 
pulse-to-pulse variability we refer to the results of Klochkov et 
al. (201 1) who have shown that for a set of pulsars the cyclotron 
line energy varies with the amplitude of individual pulses. This 
amplitude most probably reflects the instantaneous mass accre- 
tion rate. 

The luminosities for A 0535+26, 4U0115+63 and Her X- 
1 were calculated by integrating the flux of each source over a 
nearly identical energy range AL (see Table 1). For V 0332+53 
(Tsygankov et al. 2010) and GX 304-1 (Yamamoto et al. 2011) 
published values for Lx and L cyc were taken, where the lumi- 
nosity has been calculated over a slightly larger energy range 
[3-100] keV. The luminosities were computed using the source 
distances listed in Table 1 . It should be noted that uncertainties in 
the distances will create additional uncertainties in our computed 
luminosities, which have not been considered here. However, the 
typical effective uncertainties of ~ 10% in the calculated lumi- 
nosities do not strongly affect our results. 



B(R) I R 



B t \R t 



R = R t +h. 



(57) 



The corresponding variation of L cyc as a function of the emission 
height h is therefore given by 



-^cyc 

~e7 



R* + h 
R, 



-3 



(58) 



where E, is the surface value for the cyclotron energy, computed 
using Equation (l56l . Note that we have neglected the variation 
of the gravitational redshift factor, which is reasonable given the 
level of approximation employed here (Staubert et al. 2007). In 
applying Eq. (l58l to the supercritical and subcritical cases, we 
set h = h s (Eq. (1401 ) and h = h c (Eq. (IBTTi). respectively. Hence 
Eq. ( 1581 can be used to develop theoretical predictions for the 
variation of L cyc as a function of Lx for supercritical and sub- 
critical sources. 

In Fig. 3 we replot the Fig. 2 data based on a rescaling of the 
vertical and horizontal axes using the parameters L* and L cr \ t , 
respectively, which are related to each other via Eq. (153V The 
value of Lcrit used for each source is listed in Table 2, along with 
the corresponding value for the surface cyclotron energy E, ob- 
tained by comparing the theoretical variation of L cyc with the ob- 
served variation for each source. In computing L cr i t , we assume 
for all sources the same canonical neutron star mass and radius 
values M* = 1.4 M @ and R t = 10 km, and we set A = 0.1 and 
w — 1 based on the theoretical considerations discussed above. 
The vertical dashed line marked Lx = L c rit separates the sources 
into their sub- and supercritical luminosity states. It should be 
noted that, as M, and R* are also input parameters for L cr it, the 
exact positioning of the source data on the x-axis is driven also 
by the assumed canonical mass and radius values, which might 
in fact differ between the individual sources. 

Fig. 3 also includes curves representing the expected theoret- 
ical variation of the CRSF energy L cyc as a function of the lumi- 
nosity Lx, computed using Eq. (1581 with either h - h s (Eq. ( 1401 ) 
for the supercritical sources or h = h c (Eq. ( 1511 ) for the subcrit- 
ical ones. We again adopt the canonical values M* = 1.4 Mq, 
R* = 10km, A = 0.1, and w = 1, and we set the Coulomb stop- 
ping optical depth using t» = 20. The values of L» and £ are var- 
ied for each source so as to improve the agreement with the data 
(see Table 2). The parameter £ is only relevant for the supercrit- 
ical sources. The values of £ reported in Table 2 for V 0332+53 
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Table 1. Source sample characteristics. 



Source 


Instrument 


Long-term 


Pulse-pulse 


AE 


Distance 






correlation 


correlation 


(keV) 


(kpc) 


4U 01 15+63 


INTEGRAL/RXTE 


cont.dependent 1 


negative 6 


[3-60] 


8.0' 


V0332+53 


INTEGRAL/RXTE 


negative 2 


negative 6 


[3-100] 


7.5 s 


HerX-1 


RXTE 


positive 3 


positive 6 


[5-60] 


6.4 9 


A0535+26 


INTEGRAL/RXTE 


4 

no 


positive 6 


[3-50] 


2.0 10 


GX-304-1 


RXTE/Suzaku 


positive 5 




[3-100] 


2.4 11 



'Muller et al. 2011, 2 Tsygankov et al. 2010, 3 Staubert et al. 2007, 4 Caballero et al. 2007, 5 Yamamoto et al. 2011, 6 Klochkov et al. 2011, 
7 Negueruela & Okazaki 2001, 8 Negueruela et al. 1999, 'Reynolds et al. 1997, 10 Steele et al. 1998, "Parkes et al. 1980 

Notes. Our analysis is based on a combination of published data (references given) and reprocessed observational data for a number of cyclotron 
line sources [column: Source] as observed by different X-ray observatories [column: Instrument], listed here. The cyclotron line sources have been 
observed to show different types of correlation, or no correlation, of their cyclotron line energies with changes in the X-ray luminosity. 
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Figure 2. Variability of the cyclotron line energy with luminosity for different sources. The blue dashed line represents the critical 
luminosity, plotted by setting Lx = L cr j t and £ cyc = E„ where L cr ; t is evaluated using Eq. 655\ . The red dashed line represents the 
Coulomb stopping luminosity, plotted by setting Lx - L C rit and E cyc = E t , where L cou i is evaluated using Eq. d59l ). Left: observations 
on a pulse-to-pulse timescale. Right: observations on longer timescales. 



and 4U 01 15+63 are in the range ^ ~ 10" 2 - 10~ 3 , as expected 
for marginally trapped accretion columns (Becker 1998). Based 
on the results depicted in Fig. 3, we conclude that the agreement 
between the observed variation of E cyc and that predicted by the 
theoretical models developed here is reasonably close for both 
the supercritical and subcritical sources. 

A special case is 4U 01 15+63. This source was previously 
observed to show an anticorrelation (e.g. Tsygankov et al. 2007 
and references therein). The right panels of our Figs. 2 and 3 (ob- 
servations on longer timescales) include the results for E cyc and 
Lx obtained by Tsygankov et al. (2007). However, recent studies 
have shown that the presence of the anticorrelation depends on 
the choice of the continuum model (Muller et al. 201 1). 

4.2. Comparison of L x with L coul 

It is also interesting to compare the observed luminosities with 
the minimum value, L cou i, required for Coulomb stopping to de- 
celerate the flow to rest at the stellar surface, given by Eq. (1541 , 
By combining Eqs. ( f54b and ([56), we find that L cou i is related to 



E, via 

/ R t \-' 3/24 / E, y 1 ' 3 
X 1 10 km/ VlOkeV/ 

(59) 

This relation is indicated by the dashed red line in Fig. 2. 
For sources to the left of this line, we expect that the effect 
of Coulomb interactions is reduced, and the final stopping oc- 
curs via passage through a discontinuous, gas-mediated shock 
(Langer & Rappaport 1982). Hence we anticipate that the emis- 
sion region approaches the stellar surface as Lx is reduced be- 
low Lcoui. This interpretation is consistent with the observational 
data plotted in Figs. 2 and 3, which indicate that E cyc approaches 
a constant value in the limit Lx/L cou i <k 1. The left- and right- 
hand panels in Fig. 2 indicate that in fact, not much variation of 
Lcyc with decreasing luminosity Lx is observed when Lx 5 L cou i, 
as expected. 

Following our hypothesis, sources in their supercritical state 
(Lx/L c rit ~ 1) should display a negative correlation between the 



9 



P. A. Becker et al.: Bimodal spectral variation in X-ray pulsars 























Lx-Lcnt 




V 0332+53 
A 0535+26 
HerX-1 

4U 0115+63 










0.01 0.05 0.10 0.50 1.00 



5.00 10.00 



0.01 0.05 0.10 0.50 1.00 



5.00 10.00 



Figure 3. Same as Fig. 2 except with the x-axis rescaled as Lx/L c rit an d the y-axis rescaled as E cyc /E*. The values of E, and L cr j t 
used for each source are listed in Table 2. The solid curves associated with each group of data represent the theoretical predictions 



for the variation of the CRSF energy E cyc as a function of Lx for each source, computed using Eq. 
evaluated using Eq. d40l > for the supercritical sources and Eq. ( BTT i for the subcritical sources. 



with the emission height h 



luminosity and the cyclotron energy while sources in the sub- 
critical state (Lx/Z-cnt Si) should display the reverse behavior. 
V 0332+53 and Her X-l in their super- and subcritical luminos- 
ity states nicely fit into that hypothesis on both long and very 
short timescales. In particular, we note that the model param- 
eters listed in Table 2 for Her X-l are the same for both the 
long-term and pulse-to-pulse data. Hence the model developed 
here for the variation of E cyc as a function of Lx provides a ro- 
bust connection with the data across the entire range of luminos- 
ity variation for this source, which implies that the underlying 
physical mechanism of variation is the same for the long-term 
and pulse-to-pulse variations. 

The subcritical source A 0535+26 shows no significant trend 
on long timescales (Fig. 2, right), perhaps due to the fact that the 
luminosity is always close to or below the Coulomb stopping 
limit, Lcoui, which suggests that we should expect little varia- 
tion of E cyc with Lx- On pulse-to-pulse timescales, A 0535+26 
shows some suggestion of a positive correlation (Fig. 2, left), as 
expected for a subcritical source. The positive correlation sug- 
gested by the short-timescale data may reflect the fact that the 
luminosity is somewhat higher than L cou i, which places it in the 
subcritical regime according to our theory. However, we note 
that our model for the variation of E cyc with Lx does not work 
well for A 0535+26, unless we choose an unreasonably large 
value for E* relative to the observational data. We believe this 
reflects the inapplicability of our model in very low-luminosity 
sources with Lx S Lcoui- For GX 304-1, only an indication of a 
positive E cyc -Lx correlation, consistent with its subcritical state, 
can be seen, as also reported by Yamamoto et al. 201 1 . No pulse- 
to-pulse spectra are yet available for GX 304-1. 

5. Conclusions 

We have examined the hypothesis that observed bimodal vari- 
ability of the CRSF energy E cyc with luminosity Lx in accretion- 
powered X-ray pulsars reflects the dominant mode of accre- 
tion, as proposed by Staubert et al. (2007) and Klochkov et al. 



(201 1). In particular, we have derived an expression for the crit- 
ical luminosity L cr j t such that the dynamics in the supercritical 
sources is determined by the radiation pressure, and the dy- 
namics in the subcritical sources is determined by a combina- 
tion of radiation pressure and either Coulomb interactions or 
gas pressure. The detailed formula for L CI i t is given by Eq. ( [32l . 
but essentially we find that for typical neutron star parameters, 
^crit ~ 1.5 x 10 37 Z?|^ 15 erg sec -1 , where B\i is the surface mag- 
netic field strength in units of 10 12 G. 

The formula for the critical luminosity was evaluated for 5 
sources, based on the maximum value for the CRSF centroid 
energy for each source, E„, which is treated as a variable param- 
eter in our approach. The results obtained for E* are close to the 
maximum observed values for the CRSF energy. The results de- 
picted in Fig. 2 confirm that Lx > L cnt in the group 1 sources and 
Lx < Lent in the group 2 sources. The situation is less clear for 
highly variable sources with luminosity Lx that crosses over the 
line L x = L crit , such as V 0332+53 and 4U 01 15+63. These two 
sources display a negative correlation between E cyc and Lx in the 
supercritical regime, as expected, but the trend does not reverse 
as predicted by our model when Lx < L cnt . This suggests that 
these sources may actually always remain supercritical, despite 
the fact that they cross the vertical line in Fig. 3. This behavior 
can be accommodated within our model by slightly changing the 
parameters A and w in Eq. $55[ . Or, alternatively, the behavior 
of these sources could indicate that their mass and radius values 
deviate from the canonical values assumed here. 

We have developed simple physical models describing the 
quantitative variation of E cyc with Lx in the supercritical and 
subcritical sources, given by Equation (|58V with the emission 
height h set using h - h s (Equation ( 1401 )) for the supercritical 
sources (group 1), and h — h c (Equation ( BIT )) for the subcritical 
sources (group 2). In Fig. 3 the formulas we derived for E cye as a 
function of Lx were compared with the data for the supercritical 
and subcritical sources. The agreement between the theoretical 
predictions and the data suggests that our fundamental model for 
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Table 2. Theoretical parameters for each source, based on analysis of the pulse-to-pulse variability and the longer-term variability. 



Source 


Long-term 


Long-term E* 
[keV] 


Long-term L cr j t 
[10 37 erg sec" 1 ] 


Pulse-pulse 


Pulse-pulse E* 
[keV] 


Pulse-pulse L cr j t 
[10 37 erg sec -1 ] 


4U0115+63 


5.72 x 10~ 2 


17.0 


2.24 


2.14 x 10~ 2 


16.5 


2.17 


V 0332+53 


7.86 x 10~ 3 


29.7 


4.06 


1.43 x 10~ 3 


27.0 


3.67 


HerX-1 




43.5 


6.11 




43.5 


6.11 


A 0535+26 




48.0 


6.78 




48.0 


6.78 


GX 304-1 




58.0 


8.30 









Notes. The parameter <f is relevant only for the two supercritical sources, 4U 01 15+63 and V 0332+53 (see Sect. 3.1). 



the physical processes operating in these systems is essentially 
correct. 

Although the observational picture is still not complete, es- 
pecially for sources with highly variable luminosities, nonethe- 
less we believe that the emerging bimodal paradigm for the vari- 
ability of the CRSF energy with luminosity in XRBPs supports 
the hypothesis that we are seeing direct evidence for two differ- 
ent accretion regimes, depending on whether the luminosity is 
above or below the corresponding value of L cr n for the given sur- 
face magnetic field strength B*. The agreement between the the- 
oretical predictions and the observational data is rather surpris- 
ing given the level of approximation employed here. Hence we 
believe that the ideas explored here can provide a useful frame- 
work for future detailed modeling of these sources, as well as 
motivation for further observations of these and other promising 
candidates. 
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Appendix A: Effective Velocity 

The characteristic emission height in the supercritical sources 
treated in Sect. 3.1 is estimated by comparing the vertical trans- 
port time for the radiation with the mean escape time for the 
photons to diffuse out through the walls of the accretion col- 
umn. This requires a determination of the "effective velocity," 
v e ff , which is defined as the net photon transport velocity in the 
vertical direction, taking into account the competing effects of 
advection and diffusion. The former process tends to drag pho- 
tons downward toward the stellar surface, and the latter pro- 
cess tends to transport photons in the opposite direction, upward 
through the accretion column. The relationship between the ef- 
fective velocity v e g and the flow velocity v is expressed by the 
dimensionless parameter defined by 

— • (A.D 

V 

In "trapped" regions of the flow, vertical advection and diffu- 
sion are nearly balanced, and consequently v e ff <K v and £ <k 1 
(Becker 1998). Trapping tends to occur in the lower, hydrostatic 
region of the accretion column in the supercritical sources. In 
this situation, the photons tend to "hover" in a small altitude 
range until they escape through the walls of the accretion col- 
umn. Hence this effect reduces the size of the emission region in 
the supercritical sources treated in Sect. 3.1. 

The gas enters the top of the accretion column moving su- 
personically, but it must come to rest at the stellar surface. It 
follows that the flow passes through a sonic point somewhere 
in the column. The sonic point is located in the middle of the 
radiation-dominated shock, where the flow begins to decelerate 
from the incident free-fall velocity Vff (Eq. 0)). Hence the sonic 
point represents the top of the hydrostatic sinking region, where 
the radiation tends to escape, and we will therefore estimate the 
value of £ using conditions there. 
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In order to determine the flow velocity at the sonic point, it 
is useful to consider the conservation of mass and momentum in 
the hydrostatic region of the column. We have 



J = pv — const. , I = P T + pv z 



const. 



(A.2) 



where J and / denote the fluxes of mass and momentum, respec- 
tively. These two fluxes are conserved in the roughly cylindrical, 
hydrostatic portion of the accretion column. 

We can use Eqs. (lA.2b to obtain a relationship between the 
flow velocity v and the radiation Mach number, M, defined by 



M=~, 
a 



where a denotes the radiation sound speed, given by 




The result obtained is 



7 = v ( 1 + ^f 



7 

— — v c 
4 



(A.3) 



(A.4) 



(A.5) 



where v c denotes the flow velocity at the radiation sonic point, 
where M — 1 . 

The value of £ can be estimated by examining the vertical 
propagating of the photons in a radiation-dominated accretion 
column described by the exact dynamical solution obtained by 
Basko & Sunyaev (1976) and Becker (1998). This solution as- 
sumes a cylindrical geometry in the hydrostatic lower region of 
the accretion column. The total radiation energy flux in the ver- 
tical direction is given by 



E r = 4P r v + 



dP t 



n e (i\\ dz 



(A.6) 



where P r is the radiation pressure, and the first and second terms 
on the right-hand side represent advection and diffusion, respec- 
tively. We define v e ff by writing 



4P r V e « = E r . 



(A.7) 



The exact solution for the flow velocity profile v(z) in a cylin- 
drical accretion column is given by (Basko & Sunyaev 1976; 
Becker 1998) 



v(z) = v c 



14 



7 + 2e f 



1 - 



14 



7-2& 



■zJz. 



where 



3 m p c 2 

8^7 2 crj.O-|| 



(A. 12) 



(A.13) 



denotes the value of the dimensionless total energy flux (E r + 
pv 3 /2)/{Jv 2 ) at the sonic point, and 



Zst = ro 



8f t o- ± 
3cr|| 



1/2 



7 + 2e r 



In 



14 



7 - 2e c 



(A. 14) 



is the altitude of the sonic point above the stellar surface. 

Combining Eqs. ( lAlTi (IATT2t . (TATT31 . and SKJfy . we ob- 
tain after some algebra 



4v 



1 - 



7 + 2e e \ v 



14 



(A. 15) 



In particular, at the sonic point (z = z s t), we have v = v c , and 
therefore our result for £ reduces to 



2& - 1 



(A. 16) 



The dimensionless energy flux is related to the incident Mach 
number of the flow, M x , via (Becker 1998) 



- I 3 
6c ~ 2 + Ml 



(A.17) 



Combining Eqs. dA. 16b and ( IA.17I ) yields for the value of £ at 
the radiation sonic point 



1 



M 



(A. 18) 



We use this relation in Sect. 3.1, where we estimate the height of 
the emission region in the supercritical sources. 



so that v e ff represents the "effective" bulk velocity that would 
yield the correct energy flux. 

By combining Eqs. ( IA.6t and ( IA.7I ). we obtain 



Veff = v + 



c dP t 
4n e a-\\P r dz 



or, equivalently, 



p = M = i + cnip dPr 

v 4o-\\JP r dz 



(A.8) 



(A.9) 



where we have eliminated the electron number density using the 
relation J = n e m p v. The pressure P r can be expressed in terms 
of the flow velocity v by using Eqs. ( IA.2I ) to write 



(A. 10) 



p 1 (z) = I-Jv(z) 



Using this relation to substitute for the pressure P r in Eq. JA.9t 
yields 

cm p dv 

f=l . (A.ll) 

9 4a\\{I-Jv)dz 



Appendix B: Coulomb Stopping Depth 

In the subcritical sources treated in Sect. 3.2, radiation pressure 
is insufficient to decelerate the flow to rest at the stellar surface. 
In this regime, the final deceleration likely occurs via Coulomb 
interactions between the infalling plasma and the mound of 
dense gas that has built up just above the stellar surface. We can 
estimate the vertical extent of this region, and therefore obtain 
an approximation of the characteristic emission altitude in the 
subcritical sources, by computing the Thomson optical depth, 
t, measured from the stellar surface, and setting it equal to the 
value required to stop the flow, denoted by t». Nelson et al. 
(1993) carried out a detailed calculation of t„, and the final re- 
sult is presented in their equation Eq. (3.34). In this section we 
summarize the derivation. 

In a magnetized pulsar accretion column, the Coulomb stop- 
ping of the gas occurs via coupling between infalling protons 
and stationary electrons in the mound. The corresponding rate 
of change of the proton kinetic energy, 



1 2 
E P = - m p v 



(B.l) 
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is given by Eq. (3.31) from Nelson et al. (1993), which states 
that 

dE p 4-7m e e 4 , 

~~r~ — V lnA - ( B - 2 ) 

dz m e v L 

where In A c is the Coulomb logarithm and v is the velocity of 
the protons. Note that the right-hand side of Eq. (IB. 2) is posi- 
tive in our sign convention since the value of E p decreases with 
decreasing altitude. 

Using Equation ( IB. Il l to substitute for v yields the equivalent 
form 

dE 2 m„ . 

— £ = — Ann e e A In A e . (B.3) 
dz m e 

We can transform from the altitude dz to the Thomson depth dr 
using 

dr = n e cridz , (B.4) 

which yields 

dE 2 m 4^4 

~T = — In Ac. (B.5) 

Treating the Coulomb logarithm as a constant and integrating 
with respect to r, we obtain the solution 

E p (t) = E Q |l - -I j , (B.6) 
where the stopping depth, r», is defined by 

m p 47re 4 lnA c 

and the incident proton kinetic energy, Eq, is equal to the free- 
fall value, 

Eo = \ m p v\ . (B.8) 

Substituting for the Thomson cross section, <x x , in Eq. (IB. 71 ) us- 
ing 

err = ^ . (B.9) 
yields the equivalent result 

1 m„ vi 

T * = ^1 — a ~^ ■ (B.10) 

6 In Ac m e c 4 

In magnetized pulsar accretion columns, with discrete 
Landau levels, the Coulomb logarithm is given by Eq. (3.32) 
from Nelson et al. (1993), which states that 

lnA c = ln(2n max ), (B.ll) 

where the maximum excited Landau level, n max , is given by 

"max = xt; — • (B.12) 

-^*-<cyc 

Combining Eqs. ( IB. 101 ) and ( IB.1U and substituting for Vff using 
Eq. (O gives the final result, 

r , = 51.4fJ^f(-^f^— , (B.13) 
\IAM Q ) UOkmj ln(2n max ) 

in agreement with Eq. (3.34) from Nelson et al. (1993). For typ- 
ical X-ray pulsar parameters, we obtain t« ~ 20, and this is the 
value utilized in computing the characteristic emission height in 
the subcritical sources in Sect. 3.2. 
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